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Percolation clusters are random fractals whose geometrical and transport properties can be char- 
acterized with the help of probability distribution functions. Using renormalized field theory, we 
determine the asymptotic form of various of such distribution functions in the limits where certain 
scaling variables become small or large. Our study includes the pair-connection probability, the 
distributions of the fractal masses of the backbone, the red bonds and the shortest, the longest and 
the average self-avoiding walk between any two points on a cluster, as well as the distribution of 
the total resistance in the random resistor network. Our analysis draws solely on general, structural 
features of the underlying diagrammatic perturbation theory, and hence our main results are valid 
to arbitrary loop order. 



PACS numbers: 64.60.ae, 64.60.ah 



I. INTRODUCTION AND SUMMARY OF 
RESULTS 

Percolation [l[ clusters are random fractals. Due to 
their randomness, it is natural to characterize their prop- 
erties, such as geometrical and transport properties, by 
means of probability distribution functions. The perhaps 
most basic of these distribution functions is the the per- 
colation correlation function x(r — r') which measures 
the pair connection probability, i.e., the probability that 
two sites located at space coordinates r and r' are con- 
nected. For many applications, it is useful to consider 
percolation models that allow for additional processes to 
unfold on the clusters. The most prominent model of 
this kind is the random resistor network (RRN) where 
occupied bonds are viewed as resistors so that the clus- 
ters can transport electric currents. One of the typical 
questions that are addressed when studying the RRN is 
that of the total resistance R(r,r') between 2 leads lo- 
cated at points r and r' on the same cluster. Because 
the underlying cluster is random, the R(r, r') is a random 
number governed by a distribution function, P(R, r — r'). 
Experimentally, this distribution might be determined 
as follows: Measure the resistance between all pairs of 
sites with separation |r — r'| using an ohmmcter. Dis- 
card all measurements yielding an infinite value as these 
measurements probe sites belonging to different clusters. 
Prepare a histogram that collects the measured finite val- 
ues. Then, P(R, r — r')dR corresponds to the number of 
points in the histogram between R and R+dR divided by 
the number of measurements with finite outcome. A sim- 
ilar experiment has been performed, e.g., some 20 years 
ago by Rammal, Lemieux and Tremblay Q. 

Other interesting measurable quantities that are inti- 
mately related to the transport properties of percolation 
clusters include the masses (the number of bonds) of vari- 
ous fractal substructures of the clusters. The most promi- 
nent substructures of a cluster are perhaps the backbone 



(the set of bonds that carry current), the red bonds (the 
bonds that are singly connected and hence carry the full 
current) as well as the shortest, the longest and the av- 
erage self-avoiding walk (SAW) on the cluster. Their 
masses are random numbers, and the corresponding dis- 
tribution functions can be measured in much the same 
way as described above for the total resistance. 

In this paper, we study the distribution functions of 
the pair connection probability, the total resistance and 
the masses of fractal substructures by using renormalized 
field theory. Our work is based on the Harris-Lubensky 
(HL) model 0, |, [| for the RRN and its nonlinear gen- 
eralization by Harris Q. In this generalized model, the 
bonds are assumed to have a nonlinear current-voltage 
characteristic of the type V ~ I r Q. The main bene- 
fit of this generalization is that it allows to map out the 
aforementioned fractal substructures by taking the ap- 
propriate limits with respect to r: the limit r — ► — 1 + 
leads to the backbone, whereas r — * oo, r — * + and 
r — > CP lead to the red bonds, the shortest and the 
longest self-avoiding path, respectively Of course, one 
retrieves the linear RRN for r — > 1. Via these limits, 
the distribution functions of the total resistance and the 
fractal masses are all related to one single distribution 
function, namely that of the total nonlinear resistance. 
Regarding the average SAW, the story is somewhat more 
complex. Obviously, the length of the average SAW lies 
between the lengths of the shortest and the longest SAW. 
In a certain sense, the average SAW is concealed some- 
where in the discontinuity at r = 0, and it not clear how 
to extract its length by taking a limit of the nonlinear 
RRN. Recently Q, we proposed a method to calculate 
the scaling form of this length via using our real- world in- 
terpretation of the Feynman diagrams of the RRN, which 
we will explain below. In this paper, we use this method 
(which is intimately related to a correct interpretation of 
the Meir-Harris model 0) to determine the distribution 
of this length. 



2 



In their seminal paper Q, HL studied among many 
other interesting things the asymptotic form at critical- 
ity of the distribution P(R, r) of the linear total resis- 
tance for large R. Their calculation is based on the as- 
sumption that the unconditional probability P uc (R, r) for 
measuring R, which differs from P(R, r) in that infinite 
resistances are not discarded, becomes independent of r 
for large R. This assumption is problematic, however, 
because large R corresponds due to scaling to short dis- 
tances, and in critical phenomena, short distances gener- 
ically lead to extra singularities which has qualitative 
consequences [loj]. Our careful analysis, based on the so- 
called short-distance expansion as a field-theoretic tool, 
reveals that this assumption breaks down beyond mean- 
field theory, and hence the long-standing result by HL 
turns out being incorrect for dimensions d < 6. 

Before moving into the depths of field theory, we would 
like to summarize for the convenience of the reader the 
resulting picture of the various distribution functions. It 
is well established, that the pair connection probability 
obeys the scaling form 



x(r J r) = |r| 2 - d -"xHr| 1 / ,/ ) 



(1.1) 



where r is the control parameter that measures the de- 
viation from the percolation point. In bond percolation, 
e.g., t is proportional to (p c — p), where p is the bond 
occupation probability, and p c is its critical value. n and 
v are the order parameter and correlation length expo- 
nents for isotropic percolation, respectively. x(x) is a 
scaling function, whose asymptotic behavior we calculate 
for small and large x. For small x, we find 



x(x) = A Q + A x x + A 2 x 2 + A 3 x dv ~ x + 



(1.2) 



where Aq and so on are expansion coefficients. For large 
x, we obtain 



X(x) = A x (d-3+2v)"/2 exp (_ a ^) + . . . 



(1.3) 



with constants A and a. The distribution P r (R, r) of the 
total nonlinear resistance obeys at criticality, r = 0, the 
scaling form 



P r (R,T) = R- 1 fl r (R\ r \-**' v ), 



(1.4) 



where <f> r is the resistance exponent of the nonlinear 
RRN. As discussed above, the distribution function for 
the total linear resistance is implied in this result for 
r — ► 1, that of the backbone mass for r — > — 1 + and so 
on. For large s, we find that the scaling function fi r (s) 
behaves as 



n r {s) = C r ,i a-ite-W+r + c 



2 s 



-dv I <p r 



(1.5) 



where C\. r and so on are constants, and where uj is the 
so-called Wegner exponent. For small s, we get 

Q r (s) = C r fl -(«*-2+a,Waw r -") cxp (-^-"/Wr-")) 



+ 



(1.6) 



with constants C r and c r . Equations (|1.4|) to (|1.6|) also 
describe the distribution of the length of the average 
SAW; we merely have to replace 4> r by 4>q := v/vsasn, 
where vsasn is the usual scaling exponent of the aver- 
age SAW on a percolation cluster, and we also have to 
replace nonuniversal constants C r> \ and so on by other 
nonuniversal constants, say Co,i and so on. We note 
that similar asym pto tic laws were derived a long time go 
by des Cloiseaux jll| , Fisher [l2| , and McKennzie and 
Moore [l3[ for SAWs on isotropic non-disordered sub- 
strates. 

The remainder of this paper has the following outline: 
Section |TT] briefly reviews the field theoretic formulation 
of the Potts model and the HL model and its nonlinear 
generalization to provide background information and to 
establish notation. Section HTT1 treats the renormalization 
of composite fields that play a role in the short distance 
expansion. Section IIVI contains the core of our short dis- 
tance expansion. Section[V]presents the derivation of our 
final results for the distribution functions. Section I VII 
gives a few concluding remarks. 



II. FIELD-THEORETIC PERCOLATION 
MODELS BASED ON THE POTTS MODEL - A 
BRIEF REVIEW 

There are two main field-theoretic approaches for per- 
colation problems: one is based on the so-called gener- 
alized epidemic process, and the other is based on the 
g-state Potts model in the limit q — > 0. The latter is 
intimately related to the RRN and thus it is particularly 
useful for studying transport properties of percolation 
clusters. In this paper we focus on the Potts model and 
the RRN along with its nonlinear generalization. Here, 
in this section we review some aspects of the field the- 
ory of the Potts model and the RRN that are essential 
for our present work, in order to provide the reader with 
background information and to establish notation. 



A. Potts model 

1. Hamiltonian 

In the literature, the minimal model or field theoretic 
Hamiltonian for the q = (n + l)-state Potts model has 
been written in several forms. The perhaps simplest form 



G 



9k 



(2.1) 

is achieved if the (n + l)-component order parameter field 
(fik satisfies the constraint 



(2.2) 



k=0 
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because, in this case, tpk shows most directly the under- 
lying symmetry of the Potts model that arises from the 
equivalence of the n+1 states. The order parameter with 
constraint (|2.2[) constitutes an irreducible representation 
of the group of permutations of n + 1 elements, i.e., the 
symmetric group S n +i- In group theoretical language 
pu| this is the representation (n, 1) [HI. The parame- 
ter t measures the distance from the critical point in the 
phase diagram, which is given in mean field theory by 
t c = 0. Note that there is only a single expansion coef- 
ficient or coupling constant g decorating the cubic term 
since the third order direct product (n, l) 3 contains the 
identity representation only once. Hence, there is just one 
cubic invariant of the order parameter. For the model to 
describe (purely geometric) percolation, one has to take 
the limit n — > 0. 

For addressing certain questions, it can be useful to 
switch from ipk to its Fourier transform tp\ and vice versa 
via the relations 



e (fx 



+ 1 A 

1 " 



(2.3a) 
(2.3b) 



with A given by 



A = 



2irl 
n+1 



mod27r with 1 = 0, 1, ...,n. (2.4) 



The prime in Eq. (|2.3a.|) indicates that A = is excluded 
from the sum over A due to the constraint (|2.2|) . Alter- 
natively, one can admit A = if one demands that (p\ 
satisfies the constraint ^a=o = that is a consequence of 
Eq. (|2.2p . In terms of (p\, the field theoretic Hamiltonian 
for the Potts model reads 



+ f <Px<P\"P-X-X' \ ■ 

A, A' J 



'■fx 



(2.5) 



2. Renormalization 

Throughout the rest of this paper, we will draw heav- 
ily on the well established methods of renormalized field 
theory. We will outline the steps that we take as much 
as possible. Given the space-constraints of this type of 
paper, however, we have to refer to the textbooks on field 
theory [h| for general background and technical details. 

At and below the upper critical dimension d c = 6, 
mean field theory, i.e., saddlepoint-integration of func- 
tional integrals followed by naive perturbation expan- 
sions (with a physical finite cutoff), breaks down. Naive 
perturbation expansions produce IR-divcrgcncies near 
the critical point, and renormalization group methods 



must be employed instead to sum up the leading diver- 
gencies of the complete expansion. In diagrammatic per- 
turbation calculations using dimensional regularization 
instead of a finite cutoff, poles in e = 6 — d arise. These 
e-poles can be eliminated from the theory by employing 
the renormalization scheme 



vl/2 

ip -> ip = Z i <p, 

T — > f = Z~ 1 Z T T + T c 
-e„2 



9^9 = Z-WZI'*- 



u = G e fj, £ g 



(2.6a) 
(2.6b) 
(2.6c) 
(2.6d) 



where the ring ° denotes bare, unrenormalized quantities. 
u is an effective coupling constant that appears natu- 
rally in any diagrammatic calculation based on the Potts 
Hamiltonian (|2.1j) . /i is an inverse length-scale, and the 
factor a~ e makes u dimcnsionlcss in the sense that its en- 
gineering dimension vanishes. G e = T(l + e/2)/(47r) d / 2 
is a convenient numerical factor depending on space di- 
mension. The renormalization factors Z, Z T and Z u are 
known to three-loop order [lq |. As far as explicit dia- 
grammatic calculations are concerned, we will not go in 
the present paper beyond two-loop order, to which the 
renormalization factors are given by 



z = i-(„-d_ 



+ (n - 1) 



/ 5n — 11 13n — 37^\ u 



12 J 36e 



(2.7a) 



N u , . /5n-9 23n-47\ u 2 
^ = l-(»-D 7 + (n-l)(— 



(2.7b) 



. 2u 

Z u = l-(n-2) — 



7n 2 - 29n + 30 23n 2 - 99n + 1 18 \ u 2 
£ 12 J 2e' 

(2.7c) 

if one uses minimal subtraction, as we do. The critical 
value of the control parameter f c = g 4 / E S(e), where the 
Symanzik-function S(e) contains simple IR-poles at all 
values e = 2/k with k = 1,2,..., becomes formally zero 
in the e-expansion. Note that the renormalization factors 
are at any order in perturbation theory polynomials in 
n, and hence it is guaranteed that they are well behaved 
in the limit n — > 0. 



B. Nonlinear random resistor network 

In their seminal work on linear RRN, HL showed, 
based on earlier ideas by Stephen [l7| , that the field theo- 
retic Hamiltonian for the Potts model can be extended to 
describe the RRN. We refer to their field theoretic model 
as the HL model. Shortly after the HL model was intro- 
duced, Harris took up ideas by Kenkel and Straley Q 
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and generalized the HL model so that it captures nonlin- 
ear voltage-current characteristics of the type V ~ I r . 
More precisely speaking, this nonlinear RRN features 
bonds between nearest neighboring sites i and j that obey 



a generalized Ohm's law Vj — Vi 



Pi,jli,j | 



where 



Pi.j is the nonlinear bond-resistance. As is customary, we 
use here and in the following the symbol r for the non- 
linearity parameter. This parameter must not be con- 
fused with the magnitude of the space position r. The 
nonlincarity parameter always appears as a power or a 
subscript, respectively, and when read in its context, the 
meaning of r should be clear. The main benefit of the 
nonlinear RRN is that it allows to map out the important 
fractal substructures of percolation clusters by taking the 
appropriate limits with respect to r. Of course, one re- 
trieves the linear RRN for r — > 1. We will refer to the 
model introduced in Rcf. Id] as the nonlinear HL model. 



1. Hamiltonian 

To capture the RRN, one distributes the n + 1 = 
(2M+1) D states of the (n+ l)-state Potts model on a D- 
dimcnsional periodic lattice (torus) with site-coordinates 
6 = A6(h,...,j D ), where = -M, -M + 1, . . . , M. 
The combination 



MAO := 8n 



(2.8) 



is a free parameter of the theory which plays a subtle 
role. We will comment on this role as we move along. 
The order parameter field of the (nonlinear) HL model 
is given by the generalized Potts spin <pg subject to the 
constraint 



(2.9) 



It is related to its Fourier transformed counterpart via 
1 



$X : 



(2M 



i 



(2M + 1) 



l)D/2 zJ" 



ie-x, 



(2.11) 



on the form 



H 



d d r< 



J2^-x \t + wA r (\) - V< 



A, A > 



(2.13) 



where 



MA) 



D 



2x(r+l)/2 



(2.14) 



a=l 



The control parameter w is proportional to the resistance 
of the occupied bonds. The parameter r specifies the 
non-linearity of the network with r = 1 corresponding to 
the linear RRN. The HL model constitutes a Potts model 
with a specific quadratic symmetry breaking. Purely ge- 
ometrical percolation is retrieved for w — > provided one 
takes the limit D — > with M finite. In contrast, for the 
model to describe the (nonlinear) RRN, one has to take 
the limits D — > and M — > oo (keeping 9o fixed) with 
the order of the limits being of crucial importance. In 
the following, we will refer to this double limit as the 
RRN-limit. 

There are some subtle, however important, features of 
the theory associated with the RRN-limit. As long as M 
is finite, we have a model for percolation with quadratic 
symmetry breaking which can be handled by simple in- 
sertions of the operator w ip_-?A r (\)(p?. However, in 
the continuum-limit of the 0- variables, M — > oo, sums of 
the form A Ql ...A Qfc diverge like M k l 2 even if D — > 0, 
and a more careful treatment is in order 0, @. In the 
RRN-limit, the free parameter #o becomes obsolete, or in 
the language of field theory, it becomes a redundant pa- 
rameter. As such, it is possible to eliminate it via rcscal- 
ing of w. That this elimination can indeed be achieved 
hinges on the fact that in perturbation theory in conjunc- 
tion with the RRN-limit, w and A r (A) always appear in 
the combination uA r (A). This in turn implies that only 
this combination, as opposed to w or A r (A) alone, has 
a well defined engineering dimension, namely 2 like r or 
V 2 . 



where A = (Ai, . . . A^) with 



A Q — 



2irl a 



(2M + 1)A0 



mod27r/A6>, 



(2.12) 



and l a = -M, -M + 1,...,M. The prime in Eq. (f2TT0|) 
indicates that A = is excluded from the sum. In the fol- 
lowing, we will drop this prime for notational simplicity 
unless stated otherwise. Through this Fourier transfor- 
mation, the constraint (|2.9p translates into </5^_ = 0. 
Expressed in terms of (p^, the nonlinear HL model takes 



2. Connection probability and average nonlinear resistance 

The success of the HL model and its nonlinear gener- 
alization stems, at least in part, from the fact that the 
two-point correlation function 



(2.15) 



G(r,A,r, w ) = (^x( r )^A(°) 



provides for an elegant and efficient route to calcu- 
late the average resistance between two points on the 
same cluster. In the following, we will use the notation 
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G(r, t, wA r (Xj) for the two-point function, as we can be- 
cause A and w always appear in the combination wA r (\), 
see above. That the two-point function (|2 . 1 5[) is helpful 
in calculating the average resistance becomes transparent 
if one notices that it constitutes a generating function: 

G(v,r,wA r (X)) = Jx(r|C) exp(-wA r (X)R r (r\C)) 

= X(r) [exp(-wA r (X)R r (r\C)) 

J (2-16) 

where R r (r\C) is proportional to the total nonlinear re- 
sistance between an arbitrary point r and for a given 
random percolation configuration C. Here, we have ab- 
sorbed the resistance p of the conducting microscopic 
bonds into the control parameter w such that R r (r\C) oc 
p°. x(r|C) is an indicator function which is one if r is 
connected to in configuration C, and zero otherwise. 
[•••]„„ denotes the disorder average over all configura- 
tions of the diluted lattice. [•••] ov stands for disorder 
averaging conditional to the constraint that r and are 
connected. 



X (r) = [ X (v\C)] av 



(2.17) 



is the usual percolation correlation function, i.e., the 
probability for any two sites a distance |r| apart being 
connected. From Eq. Q2.16|) it follows that one can ex- 
tract the average resistance 



M r (r) = [Rr(r\C)]'„ 



(2.18) 



and its moments simply via taking the derivatives of 
G(r,T,wA r (X)) with respect to wA r (X). 

As mentioned above, the great value of the nonlinear 
RRN lies in the fact that it can be used to map out differ- 
ent fractal substructures of percolation clusters. Inspec- 
tion of the overall dissipated electric power shows readily 
that 



lim M r (r) - M bb (r) 

r — >■ — 1 + 



(2.19) 



where M&& is proportional to the average mass (number of 
bonds) of the backbone (between any two site a distance 
|r| apart, of course). Moreover, it has been shown by 
Blumenfeld and Aharony [l8| that 



where M max is the mass (length) of the longest SAW. As 
alluded to above, it is not known how to extract the mass 
(length) Msaw of the average SAW from M r . However, 
we recently proposed [U a method for calculating Msaw 
based on the RRN that does not invoke M r directly. The 
idea behind this method is exploit the real-world inter- 
pretation of the RRN Feynman diagrams, which will be 
explained in the next subsection, by placing SAWs onto 
these diagrams (for details see Ref. yj). Using this ap- 
proach, one can calculate Msaw m much the same way 
as M min and M max . 



3. Renormalization, real world interpretation of Feynman 
diagrams and scaling 

From the above, it should be clear that the renormal- 
ization factors Z, Z T and Z u for the RRN are identical to 
those for the Potts model with n + 1 = (2M+ 1)-° states. 
The control parameter w, however, needs an independent 
renormalization 



w = Z Z n 



(2.23) 



lim M r (r) ~ M red (r) 



(2.20) 



with the renormalization factor Z w having no counter- 
part in the usual Potts model. This renormalization can 
be calculated in an elegant and efficient way by using 
our real world interpretation of the HL model's Feyn- 
man diagrams [ljj [2(| Hi], [22j. This interpretation is 
based on the observation that the diagrammatic pertur- 
bation theory of the HL model in the RRN-limit can be 
formulated in such a way that the Feynman diagrams re- 
semble real RRNs. In this approach, the diagrams feature 
conducting propagators corresponding to occupied, con- 
ducting bonds and insulating propagators corresponding 
to open bonds. The conducting bonds carry replica cur- 
rents — iX. The resistance of a conducting bond is given 
by its Schwinger parameter in a parametric representa- 
tion of the propagators. Because of the interpretation of 
A (apart from the factor — i) as a replicated current, we 
refer to the space in which A takes on its values as the 
replicated current space. Likewise, we refer to the space 
in which the conjugate variable 8 lives as the replicated 
voltage space. 

When this real-world interpretation is used, the task 
of calculating Z w is by and large reduced to calculating 
the total nonlinear resistance of the conducting Feynman 
diagrams. To one-loop order, one obtains [tl |20| 



where M ror i stands for the mass of the red bonds and 



lim M r (r) ~ M min (r) 

■r^0+ 



(2.21) 



where M m ; n is the mass (length) the shortest self- 
avoiding walk (SAW). Likewise, on can show that 



lim M r (v) ~ M max (r) , 

r-»0- 



(2.22) 



Z w = 1 + yL — J dxcr(x) 

The function c r (x) is given by 

x (1 — x) 



c r (x) 



^ + {l-xf' T 



(2.24) 



(2.25) 
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Note that it vanishes for x — und x = 1 if r is positive. 
For negative r, c r (x) it takes the form 



c r {x) 



n i/k 



+ (l-x) 



l/|r 



|r| 



(2.26) 



which is equal to 1 for x = und x = 1. Here, we 
encounter another manifestation of the discontinuity at 
r = which, physically, is related to the different types 
of SAWs on the percolation clusters. 

Following standard procedures of renormalized field 
theory, one can extract from the above renormalizations a 
Gell-Mann-Low renormalization group equation (RGE) 
that governs the scaling behavior of the correlation func- 
tions. For the 2-point functions, in particular, this RGE 
leads to the well known result 

G(r, t , wA r (X)) = |r| 2 ~ d ~ r ' G(t |r| 1/ ' I/ , wA r (X) \r\^ r ^ v ) , 

(2.27) 

were G is, up to scaling factors, a universal scaling func- 
tion. The critical exponents 77 and v are known to 3-loop 
order [If}']. To 2-loop order, they are given by 



1 



206 



21" 9261 £ ' 
1 5 589 



84 37044 



r 2 . 



(2.28a) 
(2.28b) 



The resistance exponent (f> r is know for general r to 1-loop 
order up to a remaining integral Q, 
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dx c r (x) , 



(2.29) 



with c r (x) given by Eq. (|2.25|) . For specific values of r, 
(j> r is known to higher order. For the linear RRN, the 
resistance exponent <j> = (j>i is known to 2-loop order [l9l . 
HI; HI]- I n the limit r — > — 1 + relevant for the backbone, 
we calculated 4> r to 3-loop order [2(| HH [HJ . In the limit 
r — > 00 relevant for the red bonds, we showed explicitly to 
3-loop order that cf> r = 1 [20I . [2~TI . I22I as had to be expected 
from rigorous results by Coniglio [241 ]. Furthermore, we 
calculated 4> r in the limit r — > + describing the shortest 
SAW [13, E3, H3, [25| as well as the exponents governing 
the longest and the average SAW to 2-loop order 



4- Probability distribution for the nonlinear resistance - 
definition and mean field theory 

As discussed in the introduction, it is useful for a com- 
prehensive characterization of the properties of the RRN 
and the various fractal substructures of percolation clus- 
ters to go beyond just studying the average nonlinear re- 
sistance or its higher moments and to address the more 
general problem of the form of the distribution function 
of the nonlinear resistance. 

The distribution P r (R, r) of the conditional probability 
that the total nonlinear resistance between r and is R r 



given that r and are on the same cluster can be defined 
as 

P r (R,r) = [S(R-R r (r\C)}' av . (2.30) 
Comparison with Eq. (|2.16p reveals that 

1% r) := ^ F ' T '^ = r dRc~ zR P r (R, r) (2.31) 

is the Laplace transform of P r (R, r). Thus, one can cal- 
culate P r (R,r) via calculating the two point correlation 
function G(r, r, uA r (A) =: z) and the applying inverse 
Laplace transformation to IL(z, r), 



P r (R, r) 



2ni 



<7-\-i 3C 



dze zR U(z,r). 



(2.32) 



This is the main route that we will follow in the remain- 
der of this paper. 

The distributions of the various fractal masses can be 
defined by modifying Eq. (|2.30[) in an evident manner. 
These distributions can be extracted from P r (R,r) sim- 
ply by taking the appropriate limit with respect to r ex- 
cept for the distribution PsAw(£> r ) of the length L of 
the average SAW. Nevertheless, by combining the route 
described above with the method proposed in Ref. [H , we 
can calculate Psaw(L,t) in much the same way as the 
distributions of the other fractal masses. 

Before embarking on field theory, we conclude this re- 
view section by considering mean field theory At 0-loop 
order, the 2-point correlation function is given simply by 
the Fourier transform of the Gaussian propagator, 



G(r, r, z) 



exp (iq • r) 1 / t + z 

q 2 + T + z = (2n) d / 2 



x K(d-2)/2(\*\VrTz) , 



(d-2)/4 



(2.33) 



where J q is an abbreviation for ^ 2 \) d I d d Qi an( i where 
K a (x) is the modified Bessel function (Basset function) 
with index a. Using this 2-point function in conjunction 
with Eq. (|2.3ip in Eq. (|2 . 32[) and introducing the quantity 

C„W = yQ K^x)- 1 (2.34) 

which becomes unity at criticality, C Q (0) = 1, we ob- 
tain [H[ 



r((d-2)/2) 



(4i?/r 2 ) (d - 2)/2 



(2.35) 



for the resistance distribution in mean field theory. Note 
that, from this expression, one can read off the mean 
field results for the correlation length and resistance ex- 
ponents, v — 1/2 and <j) r = 1. The mean-field result 
for PsAw(-k, r ) is of the same form as Eq. (|2.35|) with R 
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FIG. 1: Form of the resistance distribution function at criti- 
cality as given in mean field theory by Eq. (|2.35fl with r = 
and d = 6. 



replaced by L. This implies the well known mean-field 
result z^saw = 1/2. 

At criticality, the distribution function takes on the 
scaling form (|1.4p . The corresponding mean field scaling 
function is plotted in Fig. Q] For large and small argu- 
ments, it behaves asymptotically as the leading term in 
Eq. (|1.6|) and Eq. (|1.5|) . respectively, with vj (<fi r — v) = l 
and du - 1 = (d - 2 + 2r))u/{2(cj) r - u)] = 2. 

As mentioned above, HL based their calculation on 
the assumption that the unconditional distribution func- 
tion P uc (R,r) = x( r )P{R, r ) becomes independent of 
r for s = R/r 2 > 1. Noting from Eq. (|2T2T)) that 
x(r) ~ | r |- d + 2 -'/ an d recalling from Eq. (|1.6p that 
fi r (s) ~ s~ 91 + • • • with gi = {du — 1)/<?V in this limit, 
we find P uc (i?, r) ~ |r| 9 with g' = 2 — rj — 1/u. In mean 
field theory, where 77 = and v = 1/2, g' vanishes so that 
P uc is indeed independent of r. Beyond mean field the- 
ory, however, g' is non-zero, and the assumption by HL 
is incorrect. Because of this shortcoming, the HL result 



HL 



(d — 2 + 7])u/(f>i for gi is not correct for d < 6, 



9i 

although it produces the correct mean-field value. 



III. REPRESENTATIONS, COMPOSITE 
OPERATORS AND THEIR RENORMALIZATION 

Now we turn to renormalized field theory to study 
the connection probability and the nonlinear resistance 
distribution beyond mean field theory. In particular, 
we will calculate their asymptotic forms in the limits 
x = rlrl 1 /^ < 1 and s = #/|r|^/ 1 ' > 1, respec- 
tively. These calculations are quite non-trivial because 
we need to determine the short-distance behavior of the 
2-point correlation function. This behavior can be an- 
alyzed systematically with help of the well established 
short-distance expansion [l(j ■ We will explain the appli- 
cation of this technique to the current problem in some 
detail in Scc. lIVI It turns out, that the composite field or 
composite operator ^(p^(r)(p_^(r) plays a pivotal role in 
the short-distance expansion. In this section we discuss 



the renormalization of this and related operators employ- 
ing group theory Our approach is guided by the 
work of Wallace and Young [l5[ , and we review parts of 
their work in a language that is adapted to the problem 
at hand. 

We are interested in the direct product of the order 
parameter with itself. As discussed below Eq. (|2.2p . the 
order parameter ipa is an object which transforms as the 
fundamental irreducible representation (n, 1) of the sym- 
metric group 5,1+1, if the invariant identity representa- 
tion (n+1) is ruled out through the constraint ^2g<p§ = 0- 
The product of the fundamental representation with it- 
self is reducible and it can be decomposed into a sum of 
irreducible representations, 

(n, 1) x ( n , 1) = (n + 1) + (n, 1) + (n - 2, 2) + (n - 1, l 2 ), 

(3.1) 

comprising the identity representation (n + 1), the funda- 
mental representation (n, 1), the symmetric tensor rep- 
resentation (n— 2, 2) and antisymmetric tensor represen- 
tation (n — 1, l 2 ). The latter is absent for a product of 
two identical fields. Applied to the product (pgLp^,, the 
decomposition leads to the identity representation 



A, 



the fundamental representation 



T y/n + 1 2 



B 



2 re VnTT 
and the symmetric tensor representation 



(3.2) 



(3.3) 



(n - l)Vn+T n(n + l) 

(3A) 

To simplify the discussion to follow, it is useful to 
switch from replicated voltage space to replicated cur- 
rent space via the Fourier transformation (|2.10j) . Then, 
the above operators take on the form 



K 



(3.5a) 
(3.5b) 



*a+a< B 



- S X+X',3> ( 3 - 5c ) 



if A and A' are non-zero. Otherwise 'fg = i? = a = 
0. Hence, there does not exist any component with zero 
current. 

As mentioned above, it will be the operator 



^a = : 2^x( r )^-X( r ) 



(3.6) 



which plays a central role in the SDE. Evidently, this 
operator satisfies A^ =0 — 0. As we will discuss in the 
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following, the renormalization of Ay is intimately related 
to that of the operators B and <I>^ y>, . Equation (|3.5p 
implies that the operators Ay, B and <fry _y are related 
for A 7^ via 



A - = — - 

A 2 A,-A 
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(3.7) 



Note that the decomposition of *4^ in irreducible parts 
contains a factor 1/n which could be potentially prob- 
lematic in the limit n — ► 0. This factor, however, cancels 
in perturbation theory. 

Form the group theoretic considerations above, we 
draw the three important conclusions, (i) The opera- 
tor tyy can mix additively with the order parameter field 
(py an d its derivatives because both transform according 
to the same irreducible representation (n, 1). (ii) The 
operators B and $^ « arc rcnormalizablc individually by 
multiplicative renormalization, because they belong to 
representations which are different from one another and 
1). (iii) The operator ^4^, which will play a central 
role in the SDE, contains both the independent renor- 
malizations of B and $^ _y and mixes these. 

Equipped with this information, we will now determine 
the renormalizations of the various composite operators. 
With the renormalization scheme (|2.6p , the renormalized 
version of the generalized Potts Hamiltonian is given by 



into renormalized correlation functions are already com- 
pletely renormalized. This implies that the renormalized 
version ^y. R °f ^a * s gi ven by 

h;R = Z u /2 h + - g (Zr-l)T ^ - - g (Z - 1) V 2 ^ 

= zwzr* £ h-j* + - Q z - 1/2 & -i)r$ x 

K 

- Z~ 1/2 -(Z -1)V 2 ^ . (3.12) 
9 

As noted above, ^ y mixes under renormalization addi- 
tively with ipy and its derivatives. 

Next, we move on to .4^ and <ty _y. Our strategy is to 
determine the renormalization of Ay with the help of that 



of Qy _y and B via relation (|3.7p . and we will switch back 
and forth between considering Ay and $>y _y to achieve 
this goal. As we know from our group theoretic consid- 
erations above, the renormalization of $^ _y is indepen- 
dent of that of B and tyy, and thus it is renormalizablc 
by multiplicative renormalization. Also because of this 
independence, we cannot simply extract the renormaliza- 
tion of $^ _y from the above. As the typical quadratic 
symmetry-breaking term in the Hamiltonian, the tensor 
$>y y, determines the so-called crossover exponent. Thus, 
we denote its renormalization factor by Z c so that 



(3.13) 



7-1/2 9 \ ^ ~ ~ ~ 
A, A' 



(3.8) 



To approach the renormalization of Ay through per- 
turbation theory, we consider its insertions into vertex 

According to the above, 



,(2) 



functions T- 

K, — K 



-p(2 

i.e., 1 .-. 



From this renormalized Hamiltonian, we can infer with- H itself is not renormalizable multiplicatively as it re- 
quires additive contributions from B, i.e., the renormal- 
ized version of Ay will be of the type 



out much effort the renormalizations of B and &y. As far 
as B is concerned, we note that if r is viewed as a function 
of r one can generate insertions of the operator —Z T B into 
renormalized correlation functions via functional differ- 
entiation with respect to r(r). Thus, the renormalized 
version Br of B is given by 



■A\-,R — % Ay 



YB 



(3.14) 



(2) 

Any result for r~ produced by diagrammatic per- 

turbation theory has to be of the form 



Br = Z r B = Z r Z 1 - ^2 V-RVR 



(3.9) 



p(2) 



\ (Sx,h + s X,-x) (! + A ^)) + B W . ( 3 - 15 ) 



Turning to ^ y, we note that the functional derivative 
with respect to </?r(r) leads to the equation of motion 



6H 



N 



5<p_Jr) 



<p(n) ■ ■ ■ 0(r N ) ) = ^(<^(ri) • • ■ ^(r,_i) 



X <p(Ti+i) • • • <p(r N )) 5(r - Ti) 



(3.10) 



where we have suppressed the current-index of most of 
the fields for notational simplicity. As a consequence of 
this equation of motion, insertions of the operator 



sn 

S( P-t 



z t t? x - zv 2 ^ + zl' 29 - J2 > ( 3 - n ) 



where A(n) and B{n) are polynomials in n which at each 
order in n have a perturbation expansion in powers of 
u containing e-poles. Summing both sides of Eq. (|3.15|) 
over A produces 

T^l^ = (I + A(n))+nB(n). (3.16) 

(2) 

We already know that is renormalized by mul- 

(2) (2) 

tiplication with Z T , i.e., Z T Y\ = _^.g is free of 
e-poles, and thus 

rl 2) _- Sfl = Z r + Z T (A{n) + nB(n)) = finite , (3.17) 
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from which Z T can be constructed order for order in u 
with polynomial factors in n. 

Equipped with this information, we now insert the 
renormalized version of A^ as given in Eq. (|3.14p into the 
renormalized two-point vertex function and determine Z' 

(2) 

and Y such that T-. s , is cured at any order in u 
from e-poles, i.e., such that 

1 H,-k;A X . r ~ k,-k;A x ' 1 1 k,-k;B 

= l fa* + (Z' + Z'A(n)) + Y 

+ Z'B(n)+Y(A(n)+nB(n)) (3.18) 

is Unite. It follows from summation over the n nonzero 
values of A that 

(Z' + nY) + (Z' + riY) (A(n) + nB{n)) = finite . (3.19) 

Comparison of Eqs. (|3.17p and (|3 . 19[) now reveals that 

Zt — . (3.20) 



IV. SHORT DISTANCE EXPANSION 



Y = 



Hence, the remaining task in determining the renormal- 
ization of A^ is to determine Z' . To this end, we sub- 
stitute the result (|3.20p into Eq. (|3. 14|) and we then note 
that 



1 



-Bp = Z' 

i 

which reveals that 



•Ax r DR 



A X 



1 



= l*x,-x 



Z' = Z c . (3.22) 

Because Z c , Z T and Y are at each order in u polynomials 
in n, it follows that in particular for n — ► that 

(3.23a) 
(3.23b) 



lim Z c = lim Z T 



lim Y = lim 

n^O n->0 



Z T — Z ( 



whereby we now have completely determined the renor- 
malization of A% in the percolation limit. Note that these 
results are valid to arbitrary order in perturbation the- 
ory. To double-check our reasoning, we also calculated 
Z c and Y explicitly in a two-loop calculation yielding 



-, (3.24a) 



Y = — 



u 




"3n-9 5n-47" 


u 2 


e 




£ 12 


4s 




"571 - 11 23?i-65" 


u 2 


+ 




e 12 


4e 



(3.24b) 



In conjunction with Eq. (|2.7b[) . the results clearly confirm 
our above general results. 

It is worth pointing out that the equality of temper- 
ature renormalization and crossover renormalization in 
the replica limit is a generic feature of all replicated field 
theories which therefore have a crossover exponent which 
is strictly equal to 1. For example, this is also true for 
the n — ► limit of the 7i-component (/j 4 -theory describ- 
ing linear polymers. Surprisingly, this important general 
feature has not been much appreciated in the literature, 
at least as far as we know. 



As discussed in Sec. Ill B 4l the distribution function for 
the nonlinear resistance can be extracted from the two- 
point correlation function G(r, r, wA r (X)), whose scaling 
form is well known to be given by Eq. (|2.27[) . For treat- 
ing the limit \r\/R^ r ^" <C 1, however, we need to know 
how the scaling function G(x,y) in Eq. (|2.27| behaves for 
small arguments. Naively, one might attempt a Taylor 
expansion of G(x, y) for small arguments which is er- 
roneous, however, because the scaling function contains 
infra-red (IR) singularities. To treat these IR singulari- 
ties properly, one has to resort to a short-distance expan- 
sion (SDE). For the </> 4 -model, this was done a long time 
ago by Brezin, Amit, and Zinn- Justin p7j . and Brezin, 
DeDominicis, and Zinn- Justin [28j]. To follow their route, 
we will briefly review some essentials of the SDE and 
its basis, the operator product expansion (OPE) [29| . 
Then we will investigate the short distance behavior of 
G(r,T,wA r (X)) by applying the SDE which will involve 
the operators . 4^ a nd B whose renormalization we dis- 
cussed in Sec. IIIII In the following, we will suppress 
the subscript R on renormalized operators for notational 
simplicity. 

The OPE is based on the fact that a product of two 
renormalized operators O a (r) and Op(r') (which may be 
plain or composite fields) at two different but close points 
can be expanded as [l(| 

O a {v)O p {v') = - r')G 7 ((r + r')/2) (4.1) 

7 

where the c Qj a. 7 (r) are pure functions of r. Of course, 
any decomposition into symmetry representations has to 
be done such that both sides of this equation remain 
consistent. The statement (|4.1[) is meaningful only when 
inserted in correlation or vertex functions and only in 
the critical domain. The same holds true for many of the 
equations in the remainder of this section, and we ask 
the reader to keep this in mind when going through this 
section. If the O a (r) are eigenoperators of the RG, i.e., 
if 



O a {lv)=r x "O a {v), 



(4.2) 



where x a is the scaling dimension of O a (r), then the 
c a p,y{r) are given by 



Car/3,7 (f) 00 



(4.3) 



The higher the scaling dimension of an operator C 7 in 
the expansion (|4.ip . the less it contributes if |r — r'| is 
small. Thus one naturally obtains an expansion for short 
distances, or SDE, if one orders the operators in the OPE 
according to their scaling dimension. 

In general, the unit operator will be present on the 
right hand side of Eq. (|4.ip . However, if one considers 
connected correlation functions, as we do, it does not 
contribute. The leading contributing operators are then 
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those which have at d c the naive dimension 4. viz. *4^(R) 
and B(R), and thus the leading contributing terms in the 
SDE of ip x (r)ip_ x (r') read 



£ x (R + r/2)v$_ x (R-r/2) 



(4.4) 



All other composite fields such as [Vcp x \/Lp_ x \ (R), 
\(pyy 2 (p _^](R), [c^^VI/ j^] (R), and their sum over A do 



at least have naive dimension 6 at d c 
true for the operator 



The same holds 



C(R) 



w 



(R) =: wC'(R) . (4.5) 



Note, that this composite operator must be renormalized 
as a whole; its rcnormalization docs not simply follow 
from the renormalization of its part [v?r</?_^] ~ At in 
the RRN-limit. This was overlooked in earlier work [15| 
and led the wrong result cf> r = 1 for the crossover expo- 
nents. However, the operator C plays a special role due 
to the double limit D — > and M — > oo. When the other 
dimension 6 operators such as ^2 X [^X^^-a]^) are m ~ 
serted into Feynman diagrams, they leave the superficial 
degree of divergence unchanged. When C(R) is inserted 
into Feynman diagrams, it reduces the superficial degree 
of divergence of any of these diagrams by 2 just as B(R) 
does. Thus, C(R) should be included on the same foot- 
ing as fi(R) as far as making insertions to avoid extra 
UV singularities if the limit |r — r'| — > is concerned. 
All other operators whose naive dimension is 6 or higher, 
can be neglected in the following. 

Having established the specific form of the SDE for 
the problem at hand, we will now turn to various RGEs 
that will eventually reveal the short distance behavior 
of the 2-point correlation function. These RGEs can be 
set up using standard procedures exploiting the fact that 
the unrenormalized theory must be independent of the 
inverse length-scale fj, introduced by renormalization [l(| , 
cf. the renormalization scheme (|2.6|) . In the RGEs for 
the current problem we encounter on various occasions 
the RG differential operator 



D M := fidf_, + (3d u + kt8 t + ( r wd u 
which contains the RG-functions 



(4.6) 



du 
0/i 



, . Slnr 
» K ( U )=M-^ — 



Cr (U) = H 



dlnw 



and the RG-function 



(4.7) 



(4.8) 



where the |o indicates that bare quantities are kept fixed 
while taking the derivatives. All these RG-functions are 



well known and explicit results for them can be straight- 
forwardly reconstructed from the renormalizations and 
Z-factors reviewed in Sec. HU 

First, let us state the RGEs for the operators A x , B 
and C, which we will need as input as we move along: 



[D^ + k] A x (R) 
[D„ + k] B(R) = 
[Dp + Cr] C'(R) = 



hyB(R) = 0, 

0, 

= 0, 



with the RG-function y(u) given by 
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65 
24 1 



(4.9a) 
(4.9b) 
(4.9c) 



(4.10) 



where we used on the right hand side our 2-loop re- 
sult (|3.24b[) for Y. Next, we set up and solve RGEs for 
the SDE coefficients c A and c B . From Eq. (|4.4[) we learn 
by using the distributivity of and by exploiting the 
RGEs glj) that 

= [X>„ + 7] (>x( R + r / 2 )^-A( R - r / 2 )) 
« [D, + 7] {c A (r)A x (R) + c B (r)B(R)) 
= A X (R) [£> M + 7 ] c A (r) + c A (r)D^A x (R) 
+ B(R) [D„ + 7] c B (r) + c s (r)£> M B(R) 
= A X (R) [D„ + 7 -K\c A (r) 
+ B(R) ([£>„ + 7 - k] cb(t) - yc A (r)) . (4.11) 

Because the operators A x and B are linearly independent , 
this implies that 

[D^+ 1 - K ]c A {r)=0, (4.12a) 
[D ll + 1 -K]c B {v)-yc A {v) =0. (4.12b) 

At the RG fixed-point it*, determined by f3(u*) = 0, these 
coupled RGEs for and c B are readily solved with the 
result 



(4.13a) 



c A (r; n) = c A (^ir) = (fj, |r|) 2 v 1/v c A {\) 

c B (r;/i) = c B (^r) = (/x |r|) 2 ~''~ lAy 

x [c B (l)-j/( U >^(l)In( M |r|)] . (4.13b) 



The additional logarithmic contribution In (/i |r|) appears 
because the operators A x and B have the same scaling 
dimension, d— 1/V, in the percolation limit n — > 0. 77 = 7* 
and v = 1/(2— «;*), where 7* = 7(u*) and so on, are the 
usual critical percolation-exponents, cf. Sec. Ill B 31 For 
completeness, and because we will use this later on, we 
also note that the resistance exponent <f> r is determined 
by the fixed point values of the RG-functions as <p r = 

"(2 -CO- 

Now we are in the position to address the short dis- 
tance behavior of G(r,T,wA r (A)). We must not naively 
simply plug the expansion (|4.4[) into the correlation func- 
tion (|2.15[) because this would generate correlation func- 
tions such as („4)*(R)) which require additional additive 
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rcnormalizations. This problem can be avoided by tak- 
ing derivatives with respect to r and w which generate 
insertions of B and C that reduce the superficial degree 
of divergence. It turns out that 3-fold differentiation is 
sufficient to suppress additional JTV-smgularities. In the 
following, we denote derivatives with respect to — r and 
— w with superscripts iV and M, e.g., 

A r (X) M G^ M Hr,T,wA r (X)) 
:= (-d T ) N (-d w ) M G(r,T,wA r (X)) 
= (^(r/2)^ x (-r/2)£(0) w C'(0) M ) , (4.14) 

where 

B(q) = f d d rB(r)e- i(iT , (4.15) 

and likewise for C. Till noted otherwise, we assume that 
N + M > 3. Then it is safe to use the SDE (JO} in 
the connected correlation function with insertions (|4.14D , 
which leads to 

A r (X) M G^ N - M Hr,r,wA r (\)) 
«(U(r) (A x (0)B(0) N C'(0) M ) 

+ c B (r) (B(p)B(0) N C' (0) m ) . (4.16) 

This equation makes it transparent that the expansion 
has the benefit of factorizing the spatial dependence from 
the dependence on the control parameters r and w. This 
feature will be important when it comes to solving the 
RGE for G (N - M \ As pointed out above, this expansion 
is valid only in the critical domain, i.e., for | 't - //j. 2 | <C 

l^ii-p 1 and \w/fJ?\ U ^ r <C |/ir| _1 . 

Next, we discuss the 2 contributions to Eq. (|4.16[) . The 
second term on its right hand side vanishes in the limit 
n — > 0. This can be seen as follows. Let us abbrevi- 
ate (B(0)B(0) N C'(0) M ) = E'x( A \(°) B ( ) Nc '( ) M ) : = 
^2'^ /(A r (A)), where we have reinstated the prime on the 
sums over replicated currents to emphasize that A = 
is excluded from these sums, c.f. the comment below 
Eq. (|2.12p . Switching to an unrestricted sum, we can 
rewrite /(A r (A)) as Ea f( A r{X)) - /(0). By express- 
ing /(A r (A)) in terms of its Laplace transform, switch- 
ing from discrete to continuous replica currents and then 
performing the integration over A with help of the saddle- 
point method, we find that J2\ fi-^r(X)) = /(0) in the 
limit D — > 0. Thus, the second term on the right hand 
side of Eq. (|4. 16|) vanishes in this limit. For discussing 
the first term on the right hand side of Eq. (|4. 1 6|) . we 
denote 

F { A N - M \ T ,wA r (X) lf i) := A r (X)- AI (A x (p)B(p) N C'{0) M ) . 

(4.17) 



An RGE for this connected correlation function is readily 
derived: 

[D M + (1 + N)k + MCr] F { a n ' m \t, wAr(X), m) = . 

(4.18) 

At the fixed point, it is straightforward to solve this RGE. 
The result is 

F a n ' m \t, wA r (X),») = +Wr 

x F A NM \t'T,l^wA r {X)4^)- 

(4.19) 

This solution needs to be augmented with a dimensional 
analysis, 

F A N ' M \T,wA r (X),n)=n d -* N+ V-™ 

xF A NM \r/^,wA r (X))/^ 2 ,l), 

(4.20) 

to account for the non-vanishing naive dimensions of its 
various ingredients. By combining Eqs. (|4. 19)) and (|4.20p 
and by exploiting freedom of choice regarding the flow 
parameter £ by setting £ = (r/fi 2 )", we obtain the scaling 
form 

F A N ' M \r,wA r (X),fx) = /-2(A r +D-2A/ 

x ( r/M y-C^)-^ 

x^(^^),(4-21) 

where F A (z) is a scaling function. For z around 
0, this scaling function is certainly analytic because the 
model is for r^0 away from its critical point, and hence 
there can be no IR-singularitics. The situation is more 
subtle for r — * because then the question arises whether 
a nonzero w makes the system noncritical or not. The 
most economical way to address this question is to cal- 
culate i^ 0,0 ' (z) in dimensional regularization to lowest 
order beyond mean-field theory. Diagrammatically, this 
calculation involves the attaching of external legs (propa- 
gators) to the usual 1-loop self-energy diagrams followed 
by an integration over the external momentum (which 
amounts essentially to a 2-loop calculation). Note that it 
is sufficient to work for N = M = in dimensional reg- 
ularization, at least formally, because this regularization 
formally sets the otherwise required additive rcnormal- 
izations equal to zero. We perform this calculation for 
the physically relevant limits of r discussed above. For 
r > 0, we find that w ^ is sufficient to make the system 
noncritical for r = 0. Thus, we expect that 

F A N ' M) (z) = K r z {du-N-i)/^-M (4 22a) 

for z — > oo, where K r is some constant. For r < it 
appears as if w ^ is not sufficient to make the system 
noncritical for r = 0. However, we encounter the problem 
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that the calculation produces a contribution proportional 
to r~ 1+£ that has no counterpart in mean- field theory, 
i.e., Fj? (z) is not rcnormalizable with the usual renor- 
malization factors for r < 0. Currently, we do not know 
how to overcome this problem or what its physical sig- 
nificance might be. Nevertheless, we expect that 



F { f' M \z) oc 1 + K r z ^-N-i)/<t> r -M 



(4.22b) 



for z — > oo if r < 0. Until the problem with the renormal- 
izability is resolved, this formula should be understood 
as a proposal, as do our final results for the distribution 
functions pertaining to r < 0. We plan to return to the 
unresolved difficulties for r < in a future publication. 

Now, we move on to the RGE for G^ N,M ^ . In its generic 
form, it reads 

[D^ + Nk + M( r ]G {N > M) (r,T,wA r (\),[i) = 0. (4.23) 

For solving this equation, it is useful to split it up into 2 
parts such that one of the parts collects the derivatives 
with respect to r and w contained in D^. At the fixed 
point u*, this procedure leads to 

+ v + Nk* + MC r ]G {N ' M) 
= k'tGW+W + C>A ? ,(A)G^ M+1 > 



CA 



k*tF 



(N+1,M) 



A* A /\\T?(N,M+t) 

( r wA r {\)F\ 



(4.24) 



where we have dropped the arguments of G^ N ' M ^ and 
so on for notational simplicity. Note that the part of the 
equation that has been moved over to the right hand side 
constitutes the leading terms in an expansion of G^ N ' M > 
in powers of r and wA r (X). Thus, after the RGE is 
written in the form of Eq. (|4.24[) . the SDE is safe for 
N + M > 2. For the right hand of the RGE (g^U, we 
can draw on the results we established above for c_4 and 
Ff' M) . We obtain 

[fid„ + r) + Nk* + MC}G {N ' M) (r, r, wA r (X),fi) 



^d-2(N+M+l) ^^2-T)-l/v (jj^ 
2 ' 



X $ 



(N.M) I wA r (X)/fI 2 



(4.25) 



where $( 7V ' M )(z) is a scaling function which sums up the 
contributions of the different F_a{z), and whose details 
are unimportant for our main quest. Any solution of the 
inhomogenous RGE consists, of course, of the general 
solution of the underlying homogenous equation and a 
particular solution of the inhomogenous equation. The 
homogenous equation is independent of r and w, and its 
solution is found to be 

G (JV ' M) (r,M)ho m = ^+^*+^*G^ M >(r,^) hom 

= d-%(N+M+\) j fid-2^)-{N+M<),r)/vQ{N,M)t£ tJlT ^ 1\ qiq 
^ ^d-2(N+l)-2M ^ < r ^-(d-2+Ti)+(N+M<f, r )/v ^ 2 g) 



As a particular solution of the inhomogenous equation, 
we find 

G( N ' M Hv, T ,wA r (\),n) pait « ^-W'+D {fl \r\f-^ 

X ( T / 2\^-l)~(N+M4>r) S (N,M) ( wA r (X)/ f l 2 \ 

\ ff i ) y (r/fx 2 )^ J ' 

(4.27) 

where S^ N ' M \z) is some scaling function determined 
by a linear first order differential equation containing 
& N > M >(z). Now, we are finally in the position to write 
down the short distance scaling form of the 2-point cor- 
relation function via adding Eqs. (|4.26j) and (|4.27j) and 
subsequent integration with respect to r and w. Intro- 
ducing the scaling variables 



x=(r/^)^\r\Y 
y=(wA r (\)/ l x 2 )(n\r\t r/u : 



(4.28a) 
(4.28b) 



we arrive at 

G(r,T,wA r (\),v) « n d ~ 2 (n |r|)~ d+2 -" {a + A h0 x 
+ A 0>1 y + A 2 , a x 2 + A ltl xy + A . 2 y 2 
+ x d »- 1 S r {y/x^)} i (4.29) 

where A$ and so on are constants, and where S r (z) is yet 
another scaling function that, in principle, can be calcu- 
lated from S^ N ' M \z). Based on Eq. (|4~22)) . we expect 
that S r (z) is analytic near z = 0, and for z 3> 1 



S r {z) cx fl(-r) + K^-^l^ 



(4.30) 



We conclude this section by asserting what additional 
contributions we get if we proceed to next order in 
the SDE [28|, i.e., if we include operators that have at 
d c a naive dimension of 6. These operators arc tA x , 

respective scaling exponents are d, d + (4> r — 1)/^, d, 
d + u> and d + 2 — where Q = (5'{u*) — e + 

■ ■ ■ is the so-called Wegner exponent. The operator 
V 2 ^ does not contribute to the SDE, because when 
inserted into correlation functions, it produces zero due 
to translational invariance: (V 2 A^(R)B(0) N C'(0) M ) = 



V 2 (A x (R)B(0) N C'(0) 



0. Therefore, the upshot 



here is that the dimension 6 operators lead to addi- 
tional analytic contributions in the brackets in Eq. (|4.29|) 
which arc of third order in x and y, such as a; 3 , x 2 y 
and so on, as well as additional non-analytic contribu- 
tions of the form x dv S^y/x^), x dv+ ^~ x S 2 {y / x^) and 
x^ d+iD >S 3 (y/x^). 



V. SHAPE OF DISTRIBUTION FUNCTIONS 

In this section, we exploit our renormalization group 
results for 2-point correlation function derived in the pre- 
vious section to determine the shape of connection prob- 
ability and the probability distribution for the nonlinear 
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resistance. In the remainder, we will suppress the inverse 
length scale \x for notational simplicity. 



A. Connection probability 



As reviewed in Sec. Ill B 21 the connection probability, 
i.e., the probability for any two sites a distance |r| apart 
being on the same cluster, is proportional to the 2-point 
correlation function of the RRN evaluated at w = 0. 
Hence, the scaling form (|2.27p of the 2-point correlation 
function implies that 



X(r ) r) = |r| 2 - d -"xHr| 1 /") 



(5.1) 



with x{x) being a scaling function that depends on the 
scaling variable x defined in Eq. (|4.28ap . In the follow- 
ing we will discuss the asymptotic forms of this scaling 
function for small and large x. 

For small x, our SDE result <|4.29[) implies that 



X(x) = A + An + A 2 x 2 + A 5 x dv ~ l + 



(5.2) 



where Aq and so on are expansion coefficients. When we 
go to next order in the SDE, we get additional contribu- 
tions proportional to x 3 , x dv und x ^ d+UJ '^ . 

Now we turn to large x, i.e., to distances |r| that 
are large compared to the correlation length £ oc | t - | ^ . 
These distances correspond to small wave vectors q, for 
which we can expand the 2-point vertex function 1^2 (q), 
i.e., the inverse of the 2-point correlation function, as 



r 2 (q,r)=r 2 (0,r){l + ?V + O(9 4 )} 



Thus, we obtain by applying the inverse Fourier trans- 
formation 



G(r,r) 



1 r e^ 1 



1 



C rf r 2 (0,r) V|r| 



4 » e"l r l/«, (5.4) 



where we have omitted higher order terms. Using 
r2(0,r) oc r 7 , where 7 = (2 — 77)1^, we get 



G(r, r) oc 



.(d-3+277)i//2 



P (d-3+2r ) ) J //2 



\d-2+ n 



-exp (-a/) , 
(5.5) 



where a is some constant. This reveals that the large- a; 
behavior of the scaling function of the connection prob- 
ability is given by 

x (x) = A a-Cd-s+auWa cxp (- ax ») + • • • , (5.6) 
with constants A and a. 



B. Resistance and fractal masses 

To establish the general scaling form of the resistance 
distribution, we recall from Sec. Ill B 41 that this distri- 
bution can be extracted from the 2-point correlation 
function via inverse Laplace transformation. Inserting 
Eq. (|2T2T)) into Eq. (|2~32"|) . we readily obtain 



P r (R,r) 



-<t>r/v 



2ni 
fT 1 B r (.R|r 



o-ioo G(X',U) v 



(5.7) 



with a scaling function E r . In the remainder, we focus 
on criticality, r = 0, where 

P r {R,r)=Rr 1 n r (R\r\-' t ' T/v ), (5.8) 

with fi r (s) = S r (s,0), and we discuss the asymptotic 
forms of the distribution when the scaling variable 

s := R\r\~^' v (5.9) 

is either small or large. 

For large s, we draw on our SDE results yielding 



G(x,y) 
G(x,0) 



= 1 + Biy + B 2 xy + B z y z 

+ x*"- 1 (S r {y/x*r) - S r {0)) + 



(5.10) 



This expansion is now inserted into Eq. (|5.7|) . In per- 
forming the inverse Laplace transformation, the analytic 
contributions in y merely produce shorts, R = 0, in the 
form of delta functions S(R) and their derivatives. What 
remains reads 



(5-3) P r (R,r) = 



2ni 



cr+'ioo 



dyx 



du-1 



x S r (y/x^) exp (yR/ \r\^'» 



( |r| 


\d-l/u ,., 




J 2m X 



a-\-ioo 



x SrirR^^/z 1 '^) exp(z) ; 



dzz (d»-l)/*r 

(5.11) 



where S r (z) 



„dv-\ 



S r {l/z' t ' r ). At the critical point, we 



obtain in particular 

n r (*) = c r ,i s-^- 1 )/^ 



C r ,2 S 



-dv/4> r 



+ C r , 3 S -(. d » + <l>r-l)/<t>r + 4 s -{d+U,)v/<t> r 

+ ■••, (5.12) 

where C r ,i and so on are constants. Here, we included the 
next higher order terms in the SDE as discussed above. 

Finally, we turn to small s. To this end, we expand 
the 2-point vertex function for small wave vectors, 

r 2 (q;r,z) 

= r 2 (0; r, z) + q 2 d q 2T 2 (q; r, z)\ q2=Q + 0(q 4 ) 

= z^f (r/z^){l + q yz^/^f 1 (r/z^)+O( q i )}, 

(5.13) 
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where we denote z = uA r (A) in conformity with our 
conventions of Sec. Ill B 41 and where /o and /i are scal- 
ing functions. At the critical point, we can replace these 
functions by constants. By dropping terms that are sub- 
dominant for Irl ^> z~ v ^ T , we obtain 



G(r,r = 0,z) = 



f *r/*r J i + flZ -2u/^ q 2 



y 



(d-3+2?j)i/ /2c/> r 



——^-ftf/*-), (5.14) 



where y = \r\/z v ^ r in accordance with Eq. (|4.28b|) . 
Feeding this intermediate result into Eq. (|5.7p . we ob- 
tain 



Pr{R,T) 



2ni 

x exp (yR/ \ r \ Vr,u - f lV 



dyy 



(d-3+2r/)iy/20 r 



(5.15) 



The remaining integration can be simplified by substitut- 
ing 

y = X*'l"v, with X<x{\T\/R'l+*f K1 - u, + r \ 

(5.16) 

which yields 

y(d-l+2rj)/2 i ra+i 
Pr ( R , r ) oc ± _ JL I dv 



R 2m 



x exp (X(«-/cW^)), 



(5.17) 



where k is some constant. For large X, we may use 
the saddle-point approximation provided that v / '<p r < 1. 
However, v j ' <p r > 1 for d = 2 and r — 1 as well as r = oo 
since v = 4/3! We finally arrive at 

O r ( S ) = C r s -(d-2+2 V ) V /2(^-u) cxp ( _ ^a-v/Ur-v)) 

+ •••, (5.18) 

where C r and c r are constants. As discussed above, it is 
straightforward to extract from these general results the 
distributions of the linear total resistance as well ass the 
distributions of the fractal masses of the red bonds, the 



backbone and the shortest and longest SAWs by taking 
the appropriate limit with respect to r. 

For Psaw{L, r), the analysisjust presented in conjunc- 
tion with the method of Ref. [8[ leads at criticality to the 
same form as given in Eqs. (|5.8[) , (|5.12p and (|5.18[) . To 
transcribe the above results to the average SAW, we just 
need to replace on the right hand side of these equations 
R by L, 4> r by </>o = v/vsaw and the constants C r .i and 
so on by other constants, say Co i and so on. 



VI. CONCLUDING REMARKS 

In summary, we calculated the asymptotic forms of the 
pair-connection probability, the distributions of the total 
resistance and fractal masses of the backbone, the red 
bonds and the shortest, the longest and the average self- 
avoiding walk between any two points on a cluster. Our 
analysis drew solely on general, structural features of the 
underlying diagrammatic perturbation theory, and hence 
our results for the form of the distributions are valid to 
arbitrary loop order, although the critical exponents fea- 
tured in these distributions are known only to finite or- 
der, of course. The distributions of the transport quan- 
tities such as the total resistance or the backbone mass 
etc. are essentially all of the same form. In this sense, 
these distributions share a high degree of universality 

As far as future directions are concerned, we think it 
is worthwhile to study the renormalizability-problem for 
r < discussed in Sec. IIVI and its implications in more 
detail. It would be interesting to have field theoretic re- 
sults for distribution functions in directed and dynamic 
percolation. Also, it would be interesting to see experi- 
mental or simulation results for distribution functions in 
any kind of percolation that might be compared to our 
results presented here or to future results. We hope that 
our work stimulates further interest in this subject. 
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